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Abstract. We present new spectroscopic and photometric observations of the young Galactic open cluster 
Westerlund f (Wd f ) that reveal a unique population of massive evolved stars. We identify ~200 cluster members 
and present spectroscopic classifications for ~25% of these. We find that all stars so classified are unambiguously 
post-Main Sequence objects, consistent with an apparent lack of an identifiable Main Sequence in our photo- 
metric data to V~20. We are able to identify rich populations of Wolf Rayet stars, OB supergiants and short 
lived transitional objects. Of these, the latter group consists of both hot (Luminous Blue Variable and extreme 
B supergiants) and cool (Yellow Hypergiant and Red Supergiant) objects - we find that half the known Galactic 
population of YHGs resides within Wdf . We obtain a mean V— My ~25 mag from the cluster Yellow Hypergiants, 
implying a Main Sequence turnoff at or below Mv = —5 (07 V or later). Based solely on the masses inferred for 
the 53 spectroscopically classified stars, we determine an absolute minimum mass of ~ f .5 x f0 3 M Q for Wd 1. 
However, considering the complete photometrically and spectroscopically selected cluster population and adopt- 
ing a Kroupa IMF we infer a likely mass for Wd I of ~ IO 5 M0, noting that inevitable source confusion and 
incompleteness are likely to render this an underestimate. As such, Wd I is the most massive compact young 
cluster yet identified in the Local Group, with a mass exceeding that of Galactic Centre clusters such as the 
Arches and Quintuplet. Indeed, the luminosity, inferred mass and compact nature of Wd I are comparable with 
those of Super Star Clusters - previously identified only in external galaxies - and is consistent with expectations 
for a Globular Cluster progenitor. 
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1. Introduction not predict what path a star of given initial mass will 

follow as it evolves from the Main Sequence through the 

Massive stars play an important role in the ecology of w -.r-r, ± * mro\ u ± t av. >/~i a-> 

K J r ° J Wolf Rayet star (WR) phase to supernova (e.g. the Conti 

galaxies, providing a major source of ionising UV radia- scenaxio . Contimm Maeder & ContiEMJ. Clearly, the 

tion, mechanical energy and chemical enrichment (Smith , , , , , , , , . , . , ,. c 

' OJ , , v most direct way to address sucn problems is to identity 

2004). However, serious gaps m our understanding of mas- , , , . , , , . ,,. un 

1 ' ' , r i an d study massive stars within clusters, since this enables 

sive stars exist throughout their hfecycle, m large part , , , , ,,. . r , „. . 

° , in. „ , us to study a coeval population at a uniform metalhcity, 

due to the rapid evolution - and hence rarity - of such , -» , a , , , A/ro 

. . , , . . where the MS turnoff mass and hence post-Mb progem- 

obiects. For example, the lack of accurate observational , , , , , , . , / , , 

J . „. . tor mass may be accurately determined (e.g. Masscy et al. 

constraints on the metalhcity dependant mass loss rates | ^QQd j) 

of such stars - particularly for short lived phases such as 

Luminous Blue Variables (LBVs) and Yellow Hypergiants Moreover, the study of clusters is important to the 

(YHGs) - restricts our ability to follow their post-Main wider galactk evolutioil) since it is generally thought that 

Sequence (MS) evolutio n. In particular, currently we can- massive gtars form preferentially ( per haps exclusively) in 

Send offprint requests to: J. S. Clark, e-mail: star clusters (e.g. de Grijs 12004J). An extreme example of 

jsc@star.ucl.ac.uk this phenomenon is observed in starburst galaxies, with 

* Based on observations collected at the European Southern star formation occurring in Super Star Clusters (SSCs), 

Observatory, La Silla, Chile (ESO 67.D-0211 & 69.D-0039) which are inferred to be several orders of magnitude more 
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massive than galactic open clusters. Such objects may 
have very different basic properties from their smaller rela- 
tives, with many authors arguing for different kinds of top- 
heavy mass distributions (e.g. M82-F; Smith & Gallagher 
2001). The relevance of such a possibility is clear when 
one considers that a large fraction of the stellar popula- 
tion of galaxies may have formed in starburst episodes in 
the distant past. If such episodes, in which very large num- 
bers of stars form, are in any way fundamentally different 
from the much more modest star formation episodes that 
we observe in the present Milky Way, our modelling of 
the history and evolution of galaxies could be profoundly 
biased. 

Unfortunately, dense clusters of massive stars are rare 
in the galaxy. The Arches & Quintuplet Galactic Centre 
clusters (with masses of ~lxl0 4 M Q ; Figer 2004) are the 
most massive known; however, observations are hampered 
by their extreme reddening. Recently, observations of the 
apparently unremarkable cluster Westerlund 1 (henceforth 
Wd 1) have demonstrated that for over 4 decades the as- 
tronomical community has possibly overlooked the pres- 
ence of a SSC in our own galaxy. 

Discovered by Westerlund ((1961), photometric 
(Borgman ct al. 119701 Lockwood 119741 and Koornneef 
\WT7\ and spectroscopic (Westerlund 113571 West87) 
surveys suggested the presence of a population of highly 
luminous early and late type supergiants. However, 
despite these observations Wd 1 languished in relative 
obscurity throughout the past decade, due in large 
part to the significant reddening (A v ~12.9 mag. Piatti 
et al. [T998) which made high resolution spectroscopic 
observations difficult. 

Motivated by the unusual radio and mid-IR properties 
of two cluster members (Clark et al. I1998II . we have un- 
dertaken an extensive program of spectroscopic and pho- 
tometric observations of Wd 1. The first results of this 
program - the discovery of an unexpected large popula- 
tion of WRs and a new LBV - have already been reported 
in Clark & Negueruela ( 2002, 2004}. Here we report on 
the remainder of the observations, specifically their use in 
constraining the properties of the population of massive 
post MS objects within Wd 1. 

The paper is therefore structured as follows. In Sect. 
2 we briefly describe the observations and the data reduc- 
tion techniques employed. Sections 3 & 4 are dedicated to 
the presentation and analysis of both spectroscopic and 
photometric datasets. Finally, we present an analysis of 
the global properties of the cluster in Sect. 5 and sum- 
marise the results in Sect. 6. 

2. Observations & Data Reduction 

Low resolution spectroscopy of cluster members over the 
red/near-IR spectral region (~6000 to 11000A) was taken 
on 2001 June 23rd to 25th from the ESO 1.52 m tele- 
scope at La Silla Observatory, Chile. The telescope was 
equipped with the Boiler & Chivens spectrograph, the 
Loral #38 camera and the #1 (night 1) and #13 (nights 



2 and 3) gratings, giving dispersions of ~ 5 A/pixel and 
~ 7 A/pixel - leading to resolutions of « 11 A and « 16A 
- respectively. 

Intermediate-resolution spectroscopy of the brighter 
cluster members was carried out on 2002 June 7th, us- 
ing the ESO Multi-Mode Instrument (EMMI) on the 3.6 
m New Technology Telescope (NTT) at La Silla. The Red 
Medium Resolution mode was used. The instrument was 
equipped with gratings #6 and #7 in the red arm and the 
new 2048 x 4096 MIT CCD mosaic (in the 2 x 2 binning 
mode). Grating #7 covered the AA6310-7835A region at 
a dispersion of ~ 0.8A/pixel, while grating #6 covered the 
AA8225-8900A range at ~0.4A/pixel. 

Data pre-reduction was carried out with MIDAS soft- 
ware, while image processing and reduction were ac- 
complished with the Starlink packages CCDPACK (Draper 
l2TiUU|) and FIGARO fShortridge lTM7|> . 

Finally photometric observations of the field contain- 
ing Wd 1 were made with the SUperb-Seeing Imager 2 
(SuSI2) direct imaging camera on the NTT; the 2 mo- 
saiced 2kx4k CCDs providing a 5.5x5.5 arcmin field of 
view. Broadband UBVRI 1 images were obtained in ser- 
vice time on 2001 August 21, with integration times for 
individual frames ranging from 1200 s (U band) to 2 s (R 
& I bands). Given the significant interstellar reddening an- 
ticipated for Wd 1 a selection of nearby red Landolt (1992) 
standards was observed. Debiasing and flat fielding (both 
dome and twilight flat field frames were obtained) were 
accomplished with the Starlink package kappa. Final 
photometry was kindly extracted for us by Peter Stetson, 
using the DAOPHOT package in IRAF 2 fStetson ri991jl . 

3. Spectroscopic classification. 

In total, 27 stars were observed at intermediate resolution 
between 6500 and 7900A and 8200 and 8900A, and an 
additional 26 stars - including the 11 WRs identified in 
Clark & Negueruela i2()02t - were observed at low resolu- 
tion between 6500 and 8900A. A finding chart for all spec- 
troscopically observed stars is presented in Fig. 1. Spectra 
are presented in Figs. 2 (new WR candidates; Sect. 3.1), 
3-7 (Super- and Hyper-giants; Sect. 3.2-3.5) and 8 (the 
unusual sgB[e] star W9; Sect. 3.3.1). 

The 6500 to 9000A window offers the prospect of 
spectral classification of moderately reddened objects for 
which traditional optical criteria are unavailable. Given 
the presence of transitions from Hi, Hei and numerous 
low - e.g. Fei, Can and Ni - and high - e.g. Hell, Cm 
and N iv - ionisation species it would appear possible 
to construct an accurate classification scheme. While au- 
thors have recognised the importance of this window (e.g. 

The U band frames produced no useable photometry of 
cluster members due the to high reddening towards Wd 1 (Sect. 
4) and so are not discussed further. 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Andrillat et al. 1995 ), comparatively few observations ex- 
ist within the literature, particularly of high luminosity 
and/or early type stars. We were therefore forced to adopt 
a twin approach for spectral classification of our spectra 
- comparison to observational data where possible (spec- 
tral type mid-B^G) supplemented by construction of a 
grid of synthetic spectra where required (OB supergiants; 
Appendix A). 

The results of the spectral classification are presented 
in Table 1, along with the original classification by West87. 
Where possible we have adopted the nomenclature of 
West87; if an individual object is resolved into two or more 
stars in our data we adopt the convention of retaining the 
designation of Wcst87, with the addition of the suffix a 
for the brightest component, b for the second brightest 
etc. We further present the J2000 co-ordinates and the 
BVRI broadband photometry for each object if available. 

Note that while we have utilised a grid of synthetic 
spectra for the classification of the OB star population 
of Wd 1 we do not claim that these represent tailored 
analyses of individual stars. Given the breadth and com- 
plexity of the spectra presented, quantitative analysis is 
clearly beyond the scope of the current work, which aims 
to provide coarse spectral classifications for the objects 
presented. 




IV44 = V N I. 
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3.1. The candidate Wolf Rayet Stars 

Clark & Negueruela |2002| reported the serendipitous dis- 
covery of 11 WR stars within Wdl in the low resolution 
2001 ESO 1.52 m dataset, comprising 6 WN and 5 WC 
stars. Analysis of the 2002 intermediate resolution NTT 
data reveals a further two candidates within the cluster. 
These are W44 and W66 (^Candidates L & M respec- 
tively, following the notation of Clark & Negueruela 2002 ) ; 
6500 to 7900A spectra for both stars are presented in Fig 
2. 

Turning first to W66, Fig. 2 reveals prominent blended 
emission in Hen 6560+Cn 6578A, Cn 6725-42+Cm 
6727-73A, Cm 7037+Hei 7065+Civ 7064A and Cn 
7236+Cm 7210-12A. Difficulty in defining a continuum 
shortwards of 7000A prevents the accurate measurement 
of line strengths for the first 3 blends, but we measure an 

EW 130±15A for the C11 7236+Cm 7210-12A blend. 

For the lowest excitation WC8 star (WR 53), the Civ 
7065A and C 11 7236A lines are of similar strength, while 
all WC9 stars have C 11 7236A >> C iv 7065A, indicating 
a classification of WC9 for W66. 

W44 demonstrates emission in the comparatively low 
ionisation lines of Ha (~ — 30A), Hei 6678 (~ —8 A) and 

7065 ( 14A), with no emission from either higher (e.g. 

Hen, Nm-iv or Ciii-iv) or lower ionisation species (e.g. 
Fell). The spectrum in the 8200 to 8900A range appears 
completely featureless (to a limit of EW ± 0.5A) and 
hence is not presented. Given the apparent magnitude of 
W44, combined with the strength of the emission lines, 
the lack of Paschen lines (emission and absorption) and 



Fig. 2. R band spectra of the 2 newly discovered WR stars 
in the cluster; Candidate L (— W44), a probable WN9 star 
and Candidate M (=W66), a sixth WCL star (WC9). 



low ionisation species, we exclude a B emission line star 
classification - such as classical Be or sgB [e] star - and in- 
stead suggest a WN identification. The lack of N iv 7116A 
precludes a classification earlier than WN9, with a clas- 
sification of WN9 rather than WN9ha more likely due to 
the comparative strength of the Hei lines. This is the first 
such star identified within Wd 1, although we note that 
the low resolution spectrum of W5 is possibly indicative 
of an even later WNL classification - this is detailed in 
Sect. 3.3.3. 

3.2. The late O & early-B supergiants 

Intermediate resolution spectra of all 13 stars in Fig. 4 
and W7, 33, 70 & 71 in Fig. 5 fulfill the criteria for OB 
spectral types defined in Appendix A and individual clas- 
sifications are listed in Table 1. The FWHM of the Paschen 
lines clearly satisfy the criterion of Caron et al. (2003) for 
supergiants, forming a continuous sequence with the later 
A-F stars in Fig. 5, which we demonstrate in Sect. 3.4 to 
also be super- or hyper-giants. 
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Table 1. Summary of the 53 cluster members with spectroscopic classifications. Column 1 provides the notation 
of West87, while Column 2 provides our additional notation if required. Columns 3 & 4 provide the Equ J2000 co- 
ordinates, while broadband BVRI photometric values are presented in Columns 5-8 where available. Columns 9 & 
10 list spectral classifications from West87 and this work, respectively. Classifications listed in italics were derived 
from low resolution spectra. In general we find an encouraging correspondence between the two works, although our 
classifications appear to be systematically earlier for those stars classified as late B-early G by West87. Importantly 
we confirm of the classification of West87's four YHG candidates, while identifying a further 2 candidates. Finally, we 
highlight one important difference; we are unable to identify any Main Sequence or giant stars, with all the OB stars 
in our study being classified as Supergiants (W6, 10, 11, 28, 33, 70 & 71 receiving revised luminosity classes). 



W87 


ID 

(new) 


RA (J2000) 


Dec (J2000) 


B 


V 


R 


I 


Spectral 
W87 


Type 
New 


2 


2a 


16 


46 


59.71 


-45 


50 


51.1 


20 


4 


16. 


.69 


14.23 


11.73 


- 


09.5 la - B0.5 la 


4 




16 


47 


01.42 


-45 


50 


37.1 


18 


.7 


14 


,47 


11.80 


9.15 


GO Ia+ 


F2 Ia+ 


5 




16 


47 


02.97 


-45 


50 


19.5 


21. 


.4 


17. 


.49 


14.98 


12.48 


- 


WNL/early B Ia+ 


6 




16 


47 


03.04 


-45 


50 


23.6 


22 


.2 


18. 


.41 


15.80 


13.16 


Bl V 


09.5 la - B0. 5 la 


7 




16 


47 


03.62 


-45 


50 


14.2 


20 





15 


,57 


12.73 


9.99 


AO la 


B5 Ia+ 


8 


8a 


16 


47 


04.79 


-45 


50 


24.9 


19 


.9 


15 


.50 


12.64 


9.89 


GO la 


F5 Ia+ 


9 




16 


47 


04.14 


-45 


50 


31.1 


21. 


.8 


17. 


,47 


14.47 


11.74 


B[e] 


sgB[e] 


10 




16 


47 


03.32 


-45 


50 


34.7 


_ 








- 


- 


B0 III 


09.5 la - B0. 5 la 


11 




16 


47 


02.23 


-45 


50 


47.0 


21 


.2 


17 


,15 


14.52 


11.91 


Bl II 


09.5 la - B0.5 la 


12 


12a 


16 


47 


02.21 


-45 


50 


58.8 


22 





16. 


,94 


13.54 


10.42 


A2 Ia+ 


A5 Ia+ 


13 




16 


47 


06.45 


-45 


50 


26.0 


21. 


.1 


17. 


,19 


14.63 


12.06 


- 


OB binary/blend? 


14 


14a 


16 


47 


05.94 


-45 


50 


23.3 


_ 








- 


_ 


- 


OB+WN blend? 


15 




16 


47 


06.63 


-45 


50 


29.7 


22 


.8 


18 


,96 


16.38 


13.75 


- 


OB binary/blend? 


16 


16a 


16 


47 


06.61 


-45 


50 


42.1 


20 


.5 


15 


,89 


12.82 


9.90 


A2 Ia+ 


A2 Ia+ 


19 




16 


47 


04.86 


-45 


50 


59.1 


22. 


.6 


18. 


,22 


15.21 


12.37 


- 


09.5 la - B0.5 la 


20 




16 


47 


04.70 


-45 


51 


23.8 


_ 








- 


- 


M6 I 


<M6 I 


23 


23a 


16 


47 


02.57 


-45 


51 


08.7 


22 


.1 


17 


,85 


14.91 


12.07 


- 


09.5 la - B0.5 la 


24 




16 


47 


02.15 


-45 


51 


12.4 


23. 





18. 


,71 


15.96 


13.24 


- 


OB binary/blend? 


26 




16 


47 


05.40 


-45 


50 


36.5 


22 


.1 


16. 


,79 


12.63 


9.19 


M2 I 


<M6 I 


28 




16 


47 


04.66 


-45 


50 


38.4 


20. 


.9 


16. 


,87 


14.26 


11.64 


B0 II 


09.5 la - B0.5 la 


29 




16 


47 


04.41 


-45 


50 


39.8 


22 


.6 


18 


,66 


16.02 


13.38 


- 


09.5 la - B0. 5 la 


30 




16 


47 


04.11 


-45 


50 


39.0 


22 


.4 


18 


,45 


15.80 


13.20 


- 


OB binary/blend? 


32 




16 


47 


03.67 


-45 


50 


43.5 


_ 








- 


- 


G5 la 


F5 Ia+ 


33 




16 


47 


04.12 


-45 


50 


48.3 


20. 


.0 


15 


,61 


12.78 


10.04 


B8 I 


B5 Ia+ 


35 




16 


47 


04.20 


-45 


50 


53.5 


22. 


.7 


18. 


,59 


16.00 


13.31 


_ 


09.5 la - B0. 5 la 


41 




16 


47 


02.70 


-45 


50 


56.9 


21. 


.3 


17. 


,87 


15.39 


12.78 




OB binary /blend? 


42 


42a 


16 


47 


03.25 


-45 


50 


52.1 


_ 








- 


- 


- 


B5 Ia + ? 


43 


43a 


16 


47 


03.54 


-45 


50 


57.3 


22. 


.8 


18 


,05 


15.22 


12.26 


- 


09.5 la - B0.5 la 


44 


L 


16 


47 


04.20 


-45 


51 


06.9 


22 


.6 


18. 


,86 


15.61 


12.52 




WN9 


55 




16 


46 


58.40 


-45 


51 


31.2 


21 


.6 


17 


,67 


15.25 


12.67 


- 


09.5 la - B0. 5 la 


56 




16 


46 


58.93 


-45 


51 


48.8 


21. 


.7 


17. 


,46 


14.81 


12.15 


- 


09.5 la - B0. 5 la 


57 


57a 


16 


47 


01.35 


-45 


51 


45.6 


20. 


.7 


16. 


,54 


13.83 


11.13 


- 


B3 la 


60 




16 


47 


04.13 


-45 


51 


52.1 


22 


.8 


18 


.50 


15.96 


13.28 


- 


09.5 la - B0.5 la 


61 


61a 


16 


47 


02.29 


-45 


51 


41.6 


21. 


.2 


17 


,16 


14.62 


12.01 


- 


09.5 la - B0.5 la 




61b 


16 


47 


02.56 


-45 


51 


41.6 


22 


.7 


18. 


,59 


16.00 


13.31 


- 


OB binary/blend? 


66 


M 


16 


47 


03.96 


-45 


51 


37.5 


_ 




19. 


,79 


16.85 


13.96 


- 


WC9 


70 




16 


47 


09.36 


-45 


50 


49.6 


21 


2 


16. 


,88 


14.10 


11.29 


B8 lab 


B3 la 


71 




16 


47 


08.44 


-45 


50 


49.3 


21 


.5 


17. 


.01 


14.06 


11.16 


B8 lab 


B3 la 


72 


A 


16 


47 


08.32 


-45 


50 


45.5 






19. 


,69 


16.59 


13.68 




WN4-5 


237 




16 


47 


03.09 


-45 


52 


18.8 


22 


.8 


17 


.49 


13.00 


9.40 


M6+ III 


<M6 I 


238 




16 


47 


04.41 


-45 


52 


27.6 


21. 


.4 


17. 


,47 


14.98 


12.45 




09.5 la - B0. 5 la 


239 


F 


16 


47 


05.21 


-45 


52 


25.0 


21 


.7 


17 


.86 


15.39 


12.90 




WC9 


241 


E 


16 


47 


06.06 


-45 


52 


08.3 
















WC9 


243 




16 


47 


07.55 


-45 


52 


28.5 














B2 la 


LBV 


265 




16 


47 


06.26 


-45 


49 


23.7 


22 


.0 


17 


.05 


13.62 


10.54 


GO Ia+ 


F5 Ia+ 




B 


16 


47 


05.35 


-45 


51 


05.0 






20 


,99 


17.50 


14.37 




WN6-8 




C 


16 


47 


04.40 


-45 


51 


03.8 
















WC8 




D 


16 


47 


06.24 


-45 


51 


26.5 
















WN6-8 




G 


16 


47 


04.02 


-45 


51 


25.2 


22 


.7 


20. 


,87 


17.75 


14.68 




WN6-8 




H 


16 


47 


03.91 


-45 


51 


19.9 


23 





18 


,55 


15.46 


12.46 




WC9 




I 


16 


47 


01.67 


-45 


51 


20.4 
















WN6-8 




J 


16 


47 


00.89 


-45 


51 


20.9 
















WN6-8 




K 


16 


47 


02.70 


-45 


50 


57.4 
















WC8-9 



Tl 



'I N f \V ™ 

^^v;^-*-*- ^fiA^. 



It-"""' 



, T .,...„ 



Fig. 3. Intermediate resolution R band (6400 to 7800 A) spectroscopy of O-F supergiant and hypergiant members of 
Wd 1. Prominent transitions of H i, He I and O I are indicated. Left panel: O and mid B supergiants; we find no evidence 
for He II absorption in the spectra. Right panel: mid-B to F supergiants in Wd 1. Weak Cn 6578/82 A absorption is 
visible in the red wing of the Ha line in the mid-B supergiants. We tentatively identify the three prominent absorption 
features between 7420-7450 A in W243 and W16a as Ni 7423.6, 7442.3 & 7468.3 A. Weak Fei lines at 7112.2 and 
6999.9 A also appear to be present in the spectrum of W12a (A5Ia+) and later objects. Possible identifications of the 
absorption features in the mid-A and later spectral types at ~7711 A and ~7748 A are Fen 7712 and Fei 7748.3 A 
respectively. 



We find W57a, 70 & 71 to have the latest spectral type 
- B3 la - of the continuous sequence of OB supergiants 
present in the sample. This classification may be assigned 
on the basis of the presence and strength of weak N I and 
O I 8446A absorption. While W7 & 33 apparently are of 
later spectral type, their spectra demonstrate several un- 
usual properties that suggest they have more in common 
with the transitional objects such as the LBV W243 and 
the YHGs; hence we defer discussion until Sect. 3.3. 

The remaining spectra display solely H I and He I pho- 
tosphcric absorption lines, implying subtypes of B1.5 la or 
earlier (T e ff > 20kK). In addition, we find no signature 
of the He II 8238A line in any object, suggesting subtypes 
no earlier than 09.5 la. Further subclassification is ham- 
pered by the poor S/N of most spectra, which precludes 
the use of the Pal6+Cin/Pal5 line ratio. 

Consequently, we may only attempt additional clas- 
sification based on the strengths of the Paschen lines. 
Measuring the strengths of the Paschen lines in the ob- 
served (and synthetic) spectra relative to W70 & 71 (B3 
la model) suggests an approximate range of spectral types 
for W61a, 60, 43a, 19, 23a, 11, 2a & 28 (earliest -> latest) 
of 09.5 la -> B0.5 la (30kK-^24kK). While the system- 
atic weakening in line strength - resulting in the absence 



of higher Paschen series lines - of stars such as W61a & 
60 compared to W2a & 28 is clearly indicative of gradu- 
ations in spectral type within this group, we refrain from 
attempting any further subclassification based on current 
data. 

Finally, while the spectra of W30, 13, 61b & 24 are 
almost completely devoid of photosphcric features and 
hence suggest spectral types of 091a or earlier, the ab- 
sence of He II 8238 A apparently precludes such a possi- 
bility. Instead we consider it likely that possible binarity 
or blending of 2 nearby objects leads to the featureless 
spectra and hence refrain from a classification other than 
generic early OB supergiant for these 4 objects. 

Although the S/N for many spectra is rather poor, 
we find no evidence for a photosphcric Ha profile as pre- 
dicted by our synthetic spectra in the temperature range 
10-34kK. A montage of the region containing the Ha line 
for those spectra with the largest S/N is presented in Fig. 
6; we find evidence for either infilling (W19, 23a, 24, 43a, 
57a 60, 61b & 70) or emission (Wll, 13, 28, 30, 61a & 
71) in all spectra. Preempting the results of Sect. 4, in the 
absence of an accurate photometric determination of bolo- 
metric luminosity, we may use the presence of Ha emission 
to infer a large luminosity and initial mass for the OB su- 
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Fig. 6. Montage of the region of the spectrum containing 
Ha for a selection of the OB supergiants. The left panel 
shows sample spectra for the 09.5 la to B3 la stars, the 
presence of either infilling or Ha emission indicating a la 
luminosity classification. The right panel shows the Ha 
profiles for the two highly luminous mid-B hypergiants 
W7 & 33, note the narrow line superimposed on a broader 
plinth of emission. 



Of the group 1 objects W44 and W13 have already 
been discussed. With an essentially featureless spectrum 
we suspect that W15 is likely a binary or blend of one or 
more OB SGs. The strong N iv emission seen in W14a is 
unexpected since this would imply an early O supergiant 
classification, at odds with the expected age of the cluster 
(Sect. 5). However, we note that W14a appears to be com- 
posed of 2 or more stars; hence an alternative explanation 
would be a blend containing at least one hitherto unde- 
tected WNE star. Clearly, future observations are required 
to resolve this apparent inconsistency. 

With a lack of O I 7774A absorption, the group 2 ob- 
jects appear to have earlier spectral types than ^B3 la. 
Indeed all those objects for which intermediate resolution 
observations are available are of spectral type 09.5 la to 
BO. 5 la - hence we assume a late O/early B supergiant 
classification for the remaining objects, with the sole ex- 
ception of W5, which, together with W7 and 33, we discuss 
in Sect. 3.3.3. 

Of the group 3 spectra W41 & 42 have only been ob- 
served at low resolution. W41 apparently shows weak O I 
absorption, although it is coupled with an unexpectedly 
weak Paschen series. Ha is clearly seen in absorption; the 
only such object in our dataset. Given the rather unusual 
nature of the spectrum we tentatively suggest the spec- 
trum may be a blend of two or more stars. By contrast 
W42a has pronounced Ha emission and appears similar 
to the mid B supergiants W7 & 33 (Sect. 3.3.3) and is 
discussed therein. 

Finally, comparison of the two epochs of observations 
of W23a and 43a are suggestive of variability, with Ha 
emission conspicuous in 2001, but apparently absent in 
2002. 



3.3. Anomalous B supergiants 

Amongst the observations of the OB supergiants within 
the cluster we find a significant number of objects with 
spectra demonstrating unusual properties. In this subsec- 
tion we detail these individually. 



pergiants. Specifically, McErlean et al. I|1999p demonstrate 
that of their sample all the most luminous (M~30-40M Q ) 
B supergiants demonstrate Ha emission. 

Fig. 7 presents a montage of the low resolution spec- 
tra of OB star candidates; where intermediate resolution 
spectra exist, these have been overplotted. Given both the 
lower resolution and S/N of the spectra, detailed classifica- 
tion is difficult. Hence we simply choose to divide the spec- 
tra into three morphological groups based on the presence 
of various photospheric features. Group 1 contains those 
stars without an obvious Paschen series and consists of 
W13, 14a, 15 & 44. Group 2 contains those stars with the 
Paschen series in absorption and consists of W2a, 5, 6, 10, 
11, 23a, 24, 29, 43a, 55, 56 & 238. Finally, group 3 consists 
of those stars with both Paschen and the O I 7774A blend 
in absorption, and consists of W7, 33, 41, 42, 70 & 71. 



3.3.1. The sgB[e] star W9 

The spectrum of W9 is presented in Fig. 8 and is dom- 
inated by H I, He I and low ionisation metallic emis- 
sion (Fen, Cai, Ni and Oi), with narrow, single peaked 
emission lines (FWHM(Ha)=125kms~ 1 ) and is notable 
for the complete lack of any photospheric features and 
the strength of certain emission features (e.g. W,\(Ha)=- 
520±30A,W A (OI)=-48±2A). To the best of our knowledge 
such an equivalent width for Ha is completely unprece- 
dented for an early type emission line star, with the ex- 
ception of the extreme systems SS433 (X-ray binary; e.g. 
Falomo et al. I1987f) and n Car (LBV; e.g. Davidson et al. 
I1999|l . Radio observations (Clark et al. 119981 Dougherty 
et al. I2()04|l reveal a strong, compact radio source to be 
associated with W9, while mid IR observations likewise 
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Fig. 7. Low resolution spectra of cluster OB supergiants (solid black lines). Where available high resolution spectra 
(from Figs. 4-6) are smoothed and overplotted (dotted red lines) on the corresponding low resolution spectra for 
comparison. Due to the low resolution an accurate determination of spectral type is impossible; hence we simply 
divide the spectra into 3 morphological classes based on the presence or otherwise of the Paschen lines and the Oi 
7774A feature; see Sect. 3.2 for further details. 



show strong emission (Fio Mm =50Jy), although the emis- 
sion mechanism (thermal dust or f-f/b-f) is still uncer- 
tain. Finally, an ISO-SWS spectrum reveals the presence 
of [O iv] emission, implying an exciting source with a tem- 
perature in excess of ^80kK must be present (Clark et al. 

imi . 

While a quantitative analysis of the combined dataset 
is beyond the scope of this work, we note that the ob- 
servational properties of W9 are unique within Wd 1, and 
apparently amongst the wider galactic population of early 
OB stars. While it fulfills the classification criteria for su- 
pergiant B[e] stars of Lamers et al. (193HJ) we suspect that 
the emission line spectrum arises at least partially in the 
compact circumstellar envelope rather than in a stellar 
wind (Clark et al. I1998[l . As such it could closely resem- 
ble the unusual object NaSt 1, which Crowther & Smith 
( 1999 ) propose to consist of a compact ejection nebula sur- 
rounding - and obscuring - an exciting WNE object. Such 
an hypothesis would naturally explain the [O iv] emission 
present in W9, although one might question the compact 
nature - and hence implied youth - of the radio nebula 
when contrasted to other LBV nebulae, given the require- 
ment for the central star to evolve through both the LBV 
and WNL phases. 

Clearly, while such an objection could be overcome by 
proposing a compact WR+LBV binary system, we also 
suggest an alternative possibility, namely that W9 is the 



result of a recent merger event. Given the severely crowded 
nature of the central regions of the cluster (Sect. 5.4) 
stellar collisions and interactions appear likely. A recent 
(~10 3 yr) interaction in which the outer H mantle of one 
or both objects were lifted off to form the compact radio 
nebula would reveal the hot interior layers of the star(s) so 
affected. This would then naturally explain the high exci- 
tation emission - until the remnant relaxed to a new equi- 
librium structure on a thermal timescale (~10 4 yrs); a sce- 
nario proposed by Clark et al. I|2000|l to explain the stellar 
and nebular morphology of the candidate LBV G25.5+0.2. 



3.3.2. The Luminous Blue Variable W243 

Clark & Negueruela (J^tHHJ present multi-epoch data that 
clearly demonstrate that W243 - initially classified as an 
early B supergiant by West87 - has evolved in the inter- 
vening decades to demonstrate a much cooler, A super- 
giant spectrum, albeit with anomalous emission features. 
The overall spectrum closely resembles that of the lumi- 
nous A supergiant IRC +10 420 (e.g. Oudmaijer 1998), 
with the notable exception of emission in the He 1 6678 
and 7065A lines. Clark & Ne gueruela ( 2004 ) advanced the 
suggestion that W243 was an LBV, with the unusual com- 
posite spectrum resulting from the formation of a pseudo- 
photosphere in a dense stellar wind. 
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Fig. 8. R and I band spectra of the sgB[e] star W9. Note that we are unable to identify the narrow emission features 
at 7134, 7233 & 7280 A. 



Subsequently, de Koter (2004, priv. comm.) suggested 
that as presented such an hypothesis would predict sig- 
nificantly stronger Ha emission (EW 100 A) than ob- 
served. However, the presence of He I emission and strong 
electron scattering wings (extending to ±10 3 kms _1 ) in the 
Ha emission profile apparently excludes a simple alterna- 
tive explanation that W243 now has the physical proper- 
ties - and in particular a rather low mass loss rate - of a 
bona fide A supergiant. 

Consequently, we suspect that a hybrid model may 
provide the correct explanation - viz. a cooling of the un- 
derlying star to a late B spectral type - still sufficient to 
produce He I emission - coupled with an increase in mass 
loss rate to simulate a later spectral type. 

3.3.3. Extreme B supergiants W5, 7, 33 & 42a 

Based on the strength of the Oi 8446A and Ni features 
present in the intermediate resolution spectra (Fig. 5), W7 
& 33 appear to be B5 supergiants. However, closer exami- 
nation shows a weak P Cygni profile in the O I 8446A line 
of both stars - to the best of our knowledge absent from 
any other mid-B supergiant, with the notable exception 
of the 2002 spectrum of the LBV W243. 

Examination of the Ha profiles also shows that they 
differ from the weak P Cygni profiles expected for B su- 
pergiants, being composed of a single narrow peak super- 
imposed on a broad emission base (Fig. 6). In this re- 
spect they again resemble the LBV W243 (Fig. 3), while 
the broad emission component is observed for the clus- 



ter member YHG W16a (Sect. 3.4). Moreover, similar 
composite Balmer line profiles are also obtained for the 
'B-supergiant LBVs' studied by Walborn & Fitzpatrick 
( 2000|. Consequently, we propose that both objects are to 
be found in a shortlived, transitional evolutionary state 
more closely related to the LBV and YHG population 
of Wd 1 than the late O/early B supergiant population. 
Indeed, both stars are a magnitude more luminous than 
any other OB supergiant within the cluster, with an ob- 
served V band magnitude directly comparable to the YHG 
population. 

Comparable objects might therefore be the highly lu- 
minous B 1.5a + hypergiants HD 80077 (Carpay et al. 
ITMT|) and HD 152236 (Van Genderen et al. HUSljl . which 
are also characterised by strong Ha emission due to mass 
loss rates an ^order of magnitude greater than la and lb 
supergiants of similar spectral type (e.g. Leitherer et al. 

Of the stars for which only low resolution spectra exist, 
we find that W42a has a comparable morphology to both 
W7 & 33, suggesting a similar classification. W5 demon- 
strates both dramatic Ha and also Cn ^7236 A emis- 
sion, previously only identified in the WCL population. 
However, the lack of Cm 9705A emision (not shown in 
Fig. 7) precludes an identification as a binary/blend con- 
taining a WCL star. Recent high resolution spectra ob- 
tained by us (Negueruela & Clark 2004, in prep) confirm 
the identification of C n emission, while also revealing the 
presence of weak He I P Cygni emission profiles. Such a 
morphology, with a lack of Niv 7116 emission, is indica- 
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tive of either a very late WN star (WN9 or later) or an 
extreme early BIa+ supergiant, the lack of a Hen transi- 
tion in the current data precluding a further classification. 

3.4. The Yellow Hypergiants 

The remaining 6 spectra in Figs. 3 & 5 are seen to be dom- 
inated by low ionisation metallic photospheric lines, sug- 
gesting spectral types later than B. Comprehensive classi- 
fication criteria and standard spectra for cool supergiants 
are provided by Munari & Tomasella ©99) . Comparison 
to their data, employing the relative strengths of the Ca II 
lines and adjacent Paschen series lines - providing the 
same diagnostic as the Call H & K versus the Balmcr 
line strength in the optical region - achieves a broad clas- 
sification for B8-F8 supergiants. This may then be further 
refined by consideration of the presence and strength of 
the Ni transitions between 8650 to 8750A for A super- 
giants and the appearance of Fei and Sil lines for the 
F supergiants in our sample. When combined with the 
(lower resolution) observations of Cenarro et al. (2001) 
we estimate a classification accuracy of ±2 spectral sub- 
types for our current observations, from which we identify 
two A and four F stars (Tables 1 & 2). 

The narrow line widths, forming a smooth progres- 
sion from the OB supergiants, together with the pro- 
nounced Ni absorption features (c.f. Figs. 17-25 of Munari 
& Tomasella 1999) clearly justifies at least a supergiant 
classification for the 6 stars in question. However, deter- 
mination of the bolometric luminosity of A-G stars may be 
made directly via the strength of the O I 7774A blend; a re- 
lationship first identified by Merrill ( 1925$. We make use of 
the most recent calibration of this relationship by Arellano 
Fcrro et al. (2003; their Eqn. 2); which is calibrated from 
Mv=+0.35 to —9.5 mag. The results are presented in 
Tabic 2, with four stars (W4, 12a, 16a & 265) apparently 
demonstrating absolute visual magnitudes in excess of the 
current calibration for the relationship. Consequently, we 
simply limit ourselves to the conservative statement that 
these stars have My ~ —9.5 (log(L*/L Q )~5.7); placing 
them at the empirical Humphreys-Davidson limit for cool 
hypergiants (under the assumption of negligible bolomet- 
ric correction). We further note that both the presence of 
chemically processed material at the stellar surface and/or 
infilling of the blend due to emission from a stellar wind 
will only cause a systematic underestimate of the true 
luminosity; such an effect is apparently observed for the 
YHG p Cas. Therefore, we conclude the six A-F stars in 
Wd 1 have luminosities that are at least directly compa- 
rable to those of the field populations of YHGs found in 
both the galaxy and the Magellanic Clouds. 

However, as described by de Jager Ijl998(l a high lumi- 
nosity (My > -7) is not sufficient to designate a star as a 
bona fide YHG. Rather, he defines the 'Keenan-Smolinski' 
criteria for the spectroscopic classification of hypergiants; 
(i) the presence of one or more broad components of Ha 
and (ii) absorption lines that are significantly broader than 



Table 2. Equivalent widths for the Oi -7774A blend for 
the six yellow hypergiants in Wd 1, with absolute visual 
magnitude and resultant luminosity determined directly 
from the relationship of Arellano Ferro et al. ( 2003 , formal 
error of ±0.4 mag.). 



YHG 


Spec. 


EW(Oi) 


M„ 




Type 


(A 




4 


F2Ia+ 


3.0 


-10.1 


8 


F5Ia+ 


2.1 


-8.7 


12a 


A5Ia+ 


2.8 


-9.9 


16a 


A2Ia+ 


2.7 


-9.8 


32 


F5Ia+ 


2.3 


-9.1 


265 


F5Ia+ 


2.6 


-9.7 



those of la stars of similar spectral type and luminosity. 
Both criteria are designed to identify those stars with an 
enhanced mass loss rate, which de Jager (1998) introduces 
as the defining criteria of YHGs. Comparison of our spec- 
tra to those of bona fide YHGs demonstrates that their 
line widths are indeed consistent with such a classification 
(H. Nieuwenhuijzen priv. comm. 2003). 

Between 2001-02, Ha has only been observed in emis- 
sion in W16a, although West87 report emission in both 
W16a and W265, suggesting that it may be a transient 
phenomenon. Indeed, transitions from absorption to emis- 
sion have been observed in p Cas (Lobel et al. 120031 . sug- 
gesting that the lack of Ha emission in the majority of 
the stars considered here may not disqualify them from 
a hypergiant classification. Moreover, we note that radio 
observations of the cluster (Dougherty et al. 200U reveal 
all six to have radio counterparts which we suggest to be 
the result of the ambient UV radiation field provided by 
the hot star population of the cluster ionising their stellar 
winds (given that the stars themselves are likely too cool 
to ionise their own winds) . We interpret this observation as 
an alternative diagnostic for mass loss from the stars and 
hence propose a hypergiant classification, noting that all 
six objects lie on the border or within the 'yellow void', a 
region of the H-R diagram occupied by known YHGs (e.g. 
Fig. 1 of dc Jager irU9"5|l . 

3.5. The Red Supergiants; W20, 26 & 237 

Finally we turn to the three M supergiants within Wd 1. 
Unfortunately, our current low resolution data do not 
allow us to improve on the classification of West87 for 
these objects. However, unlike West87, who concludes that 
W237 is seen in chance projection against Wdl, we sug- 
gest that it is in fact a bona fide member of the cluster. We 
base this conclusion not only on the comparable optical- 
mid- IR fluxes of the three M stars, but, more compellingly, 
the fact that all three are strong, spatially resolved radio 
sources (Dougherty et al.[2004). To the best of our knowl- 
edge no other RSGs have similar radio properties (e.g. 
Clark et al. I1998[) : hence we consider it likely that this 
phenomenon results from membership of Wd 1, in a man- 
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ner identical to that described above for the YHGs (Sect. 
3.4). 

The spectrum of W26 also demonstrates a strong dou- 
ble peaked emission line at 6564/82A (EW=-74A) that 
we attribute to a blend of Ha and [Nil] 6583. 5A emis- 
sion, noting that West87 also identifies Ha in the spectrum 
of W26 and W20. Two narrow emission lines at 9062 & 
9525A are also observed in the spectrum of W26. Owing to 
the difficulty of defining continua in these regions of the 
spectrum we have been unable to determine an EW for 
either line, but tentatively identify them as [S in] 9069.4 
& 9532. 5A. W26 is known to be identified with extended 
radio and mid-IR emission (Clark et al. I1998f) with an ap- 
parent bow shock morphology and we consider it likely 
that the emission arises in this nebular material. Identical 
emission lines are also observed towards two stars ~1.5 ar- 
cmin to the West (with colours consistent with member- 
ship of Wd 1), which are coincident with a an extended 
region of radio emission (Dougherty et al. I2004fl . Thus 
we suggest a 'nebular' rather than 'stellar' origin for the 
emission lines, although we note that the nebular material 
associated with W26 likely results from the interaction of 
its stellar wind with the intracluster medium/wind. 

4. Photometric results 

In Fig. 9 we present V vs. (V— R) and (V— I) colour mag- 
nitude plots for the 1044 stars within the 5.5x5.5 ar- 
cmin field containing Wd 1, with photometric errors de- 
rived from daophot presented for BVRI bands in Fig. 
10. Several distinct populations of stars are clearly de- 
lineated in both colour magnitude plots. Through com- 
parison to the similar plot for the highly reddened clus- 
ter Pismis 23 (Piatti & Claria 2 0U2|) we identify a long 
blue tilted Main Sequence (MS) for field stars from V=12- 
21 mag. and (V— R)=0.6-1.4 and a vertical plume centred 
on (V— R)=1.2, which we identify with a corresponding 
population of Horizontal Branch stars. A second, more 
sparsely populated MS is apparent at (V— R)~1.2-1.6 
- suggesting the presence of an intervening absorber 3 . 
Finally, we find a large population of highly reddened 
stars with (V— R)>2 ((V— 1)>4) that we propose as clus- 
ter members. 

Photometry for those cluster members with spectro- 
scopic identifications is provided in Table 1 . Unfortunately 
photometry of 12 spectroscopically identified cluster mem- 
bers was unavailable due to severe blending or the star 
falling in the gap between the two detectors. The LBV 
W243 is one particularly regrettable example of the latter 
failing, as it had not been recognised as a variable star 
when the observations were made. In total we have cur- 
rently identified 201 probable cluster members via a com- 
bination of spectroscopic and photometric observations. 

3 The photometric data of West87 - which were obtained 
for a greater field of view - reveal a much richer second MS, 
suggesting the obscuring cloud is of limited spatial extent, ap- 
parently covering a region not much larger than the cluster. 



4.1. The Reddening & Distance to Wd 1 

As Piatti et al. l|TM5|l . we find a scatter in the optical 
colours of cluster members that clearly exceeds the pho- 
tometric errors, and implies significant differential extinc- 
tion across the cluster. However, we note one essential dif- 
ference between Piatti et al. and our results - we find no 
evidence for a Main Sequence in our current photometric 
data set. We may use our spectroscopically identified OB 
stars to study the reddening to Wd 1 - given that the op- 
tical colours of all early OB supergiants are rather similar 

- finding no systematic gradients in reddening across the 
cluster. In Fig. 11 we plot the (V— R) vrs. (V— I) colours 
of the early and mid-OB supergiants to demonstrate the 
extent of the differential reddening. We find E(V-R) and 
E(V-I) to range from ~2.6-3.2 and ~5.3-6.2 with median 
values of 2.76 and 5.65 respectively 4 . We also plot least 
square fits to the data, allowing us to derive the redden- 
ing vector for the cluster. The resultant reddening vector 
does not pass through the locus for unreddened OB super- 
giants, indicating a non standard extinction law towards 
Wd 1; a fact we return to below. 

Given that a spectroscopic discriminant exists for the 
intrinsic luminosity of the YHGs, we may use these stars 
to determine V-My to the cluster. Again, due to differen- 
tial reddening we find this to range from 24-26.5 mag. - the 
latter value for the isolated star W265 to the North of the 
cluster (noting a formal error of 0.4 mag. due to the un- 
certainties in the calibration of the My - EW0/7774 rela- 
tionship). The mean and median values of V-My are both 
^25.3 mag., which we adopt for the remainder of the pa- 
per. Piatti et al. fiMSj l derive a value of V-Af v =23.8±0.3 

- at the lower range of values we find - via isochrone fitting 
to the sequence of OB supergiants within Wd 1, which 
they mistook for a bona fide MS. They noted the diffi- 
culty of such a task given that their erroneous MS is both 
^vertical and significantly broadened in Wd 1 due to dif- 
ferential reddening. Their vertical placement of the 4 Myr 
isochrone was therefore accomplished via reference to the 
presence of the as then currently identified supergiants, 
for which Piatti et al. (1998) significantly underestimated 
their absolute visual magnitudes, leading to a correspond- 
ing underestimate of V— My. 

Two other methods employed by Piatti et al. (1998) 
to infer V— My - and hence the reddening and distance 

- for Wd 1 were comparison of their photometry to the 
values of My quoted by West87 for the brightest clus- 
ter members and a comparison to template clusters. The 
former suffered from the significant underestimation - by 
>1 mag. on average - of the relevant Mys by West87, while 
the latter invokes a mistaken comparison to the 10-15Myr 



4 Note that since the differential reddening is a real physical 
effect we present a range of values rather than a formal error 
derived from the standard deviation of the datatset and we 
choose to use the median rather than the mean value of colour 
excess since we find that latter to be skewed by a few stars 
with particularly high reddenings. 
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Fig. 9. Plot of V band magnitudes versus the (V — R) and (V — I) colour indices for the 5.5x5.5 arcmin field centred 
on Wd 1. Objects with spectroscopic classifications are indicated with filled symbols with red, yellow and blue squares 
for RSGs, YHGs and OB SGs respectively (the sgB[e] star Wd 9 is indicated by a circle). WRs are indicated by 
magenta symbols, with squares WC and triangles WN stars; the 2 candidate WNL stars W5 & 44 are included in this 
category. Note that photometry is lacking for a number of cluster members such as the LBV W243, YHG W32, the 
RSG W20 and 6 WRs, while the two anomalously luminous OB SGs are the extreme B5 Ia + objects W7 & 33. 



template open clusters NGC 457 and NGC 884; both cases 
leading to underestimates of V— My. 

Adopting V— My ~25.3mag. to Wd 1 clearly demon- 
strates why no MS objects have yet been spectroscopically 
identified. An 07 V star (Sect. 5.1) with an absolute mag- 
nitude of V=— 4.7 would have an apparent magnitude of 
V~20.6 mag, at the limits of detectability for our current 
photometry and beyond the reach of our current spec- 
troscopy. Consequently, we are led to the conclusion that 
the majority of cluster members currently identified are 
post-Main Sequence objects. 

Finally, we may attempt to determine the distance to 
the cluster by individually dereddening the YHGs, for 
which we have a spectroscopic luminosity discriminant 
(which is lacking for the OB supergiants). This yields a 
mean distance of ~5.5 kpc and a mean A„ — '11.6 mag. 
We regard this distance estimate as an upper limit for 
three reasons. Firstly increasing the distance would yield 
correspondingly larger luminosities for both the YHGs 
and RSGs, placing them above the empirical Humphreys- 
Davidson (HD) limit for cool stars. Indeed such a large 
distance - possibly biased by an overestimate of the lumi- 
nosity of the isolated YHG W265 - already places W16 
~0.5 magnitudes above the HD limit. Secondly, the dis- 
tance determination assumes a standard reddening law 
which we know to be incorrect. From the median E(V- 



I)~5.65 we find A„=11.0, via the ratio E(V-I)=1.6E(B- 
V) given by e.g. Fitzpatrick l(T9^9) and A 1 ,=3.1E(B-V); 
appropriate for a standard reddening law. However, from 
the OB supergiants we measure a median E(B-V)~4.35, 
resulting in an At, = 13.6. Therefore we conclude that the 
above methodology results in an underestimate of the true 
reddening, leading in turn to an overestimate of the dis- 
tance to Wd 1. The final reason to favour a distance of 
<5.5 kpc is that anything greater would place Wd 1 within 
4 kpc of the Galactic Centre, a region underpopulated by 
Hn regions from which it might have formed (e.g. Fig. 5 
of RusseilEESI). 

A lower limit to the distance may be determined from 
the lack of radio detections of the Wolf Rayet population 
of Wd 1 (Dougherty et al. in prep.). Assuming a mass loss 
rate of ~lO _5 M0yr _1 and a wind velocity of 3000kms _1 
for the WRs, we would have obtained 3<r detections even 
at a distance of 2kpc (Dougherty priv. comm. 2004). 
Consequently we adopt 2kpc as a lower limit to the dis- 
tance of Wd 1 5 . 

5 We note that even this conservative estimate for the min- 
imum distance to Wd 1 is significantly in excess of the 
1.1±0.4 kpc finally adopted by Piatti et al. JEHU- Such a 
low distance estimate resulted from the errors in determining 
V— Mv described above and also the value adopted for the 
E(B-V)/E(V-I) ratio. Specifically, Piatti et al. iflM^I employ 
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4.2. The structure of Wd 1; spatial extent and mass 
segregation 

In Fig. 12 we present a plot of the positions of poten- 
tial cluster members based on the colour cut offs above. 
Previous analysis by Piatti et al. (19(21) na cl suggested 
that Wd 1 has a core+halo structure with radii of ^36" 
and ^72" respectively, with a radial light profile that di- 
verges from that expected for a King profile. From Fig. 
12 we find that ~50% of the cluster members are found 
within a ^50 arcsec diameter 'core', located approxi- 
mately at R.A.=16 47 04 5=-44 09 00 (J2000), corre- 
sponding to <1.2pc for a distance of <5.5 kpc. An appar- 
ent deviation from spherical symmetry is present, with an 
overabundance of stars to the South West - however due 
to the requirement to site the cluster away from the gap 
between the CCDs we lack information on the extent of 
the cluster in this direction. To the North East we find a 
few candidate cluster members out to radii of >3arcmin 
(as measured from our cluster 'core'). 

However, these values were determined solely from 
analysis of the distribution of the cluster supergiants, thus 
excluding both the more massive and evolved WR compo- 
nent, and also the currently undetectable Main Sequence 
(M<30M Q ; Sect 5.1). Consequently, we refrain from draw- 
ing firm conclusions on the extent and cluster density pro- 
file. We do not expect the positions of the currently ob- 
served post-MS population to accurately represent the un- 
derlying stellar population, as many young clusters such 
as R136/NGC 2070 show evidence of strong mass segrega- 
tion (e.g. Meylan 1993, Schweizer I2004f) . Indeed, by anal- 
ogy with other young massive clusters we might expect the 
density profile of Wd 1 to be described by a EFF (Elson et 
al. HWjl King-like profile ( Schweizer IStMl Larsen I2UU3|) ; 
this will be addressed in a future paper in which we anal- 
yse VLT (NAOS CONICA & FORS1) imaging data. 

5. Discussion 

Spectroscopic and photometric analysis of the upper 
reaches of Wd 1 reveal a unprecedented rich population 
of massive post-MS objects. Prior to this study, the three 
Galactic Centre clusters had hosted the richest popula- 
tion of such stars in the galaxy. In Table 3 we summarise 
the spectroscopic findings for Wdl and the populations 
of the Quintuplet, Arches and Galactic Center clusters. 
Clearly, Wd 1 rivals, if not exceeds the stellar content of 
each of these clusters. While OB supergiants and WRs are 
found in all four clusters, a rich population of cool tran- 
sitional objects is unique to Wdl; indeed, the population 
of YHGs in Wdl is equal to that found for the Galaxy or 
the LMC. Given the apparent brevity of this evolutionary 

the calibration of Dean et al. {1978), while we have adopted 
the more recent calibrations of Rieke & Lebofsky 1198511 and 
Savage & Mathis (1979), which have subsequently been con- 
firmed by other studies (e.g. Cardelli et al.[l989 and Fitzpatrick 
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Fig. 10. Magnitude errors from DAOPHOT for the cluster 
field. 




V-R 

Fig. 11. Plot of the (V-R) and (V-I) colours for the 09.5 
to BO. 5 (filled symbols) and B3 to B5 supergiants (open), 
demonstrating the extent of the differential reddening 
across the cluster. Best fits to the reddening vector of 
the early and early+mid B supergiant datasets are also 
plotted and are found to be almost indistinguishable. 

phase, inferred from the rarity of these stars, such a popu- 
lation suggests that we are still incomplete for both their 
likely progenitors - OB supergiants - and descendants - 
WR stars. 

Comparison to the current photometric dataset sup- 
ports this conclusion. Assuming an average M\ < —7 
and V—My ~25.3 mag. yields a total of 23 candidate 
OB supergiants with V<18 mag. (excluding known cool 
transitional objects); we have spectroscopically confirmed 
the identity of only ~60% of this sample. Moreover, we 
note that a further 12 of the spectroscopically confirmed 



16 



J. S. Clark et al.: The stellar content of Wd 1 



Ti- 


1 , , 1 , I 1 


- 

!•' v. ■ ■ ■ « • 

IT ■ ■ ■ ■ , . 


en 


. : . ft. 

• ' * . 




lj m 


■ + . 

m . 


ft j * 

• V •%« ' ■ • 

• .... ■ ■■•■* . ■ ■ 

I * ".-»'■ J 


Q 

LO 

lo 

"rO 


: '.«>-* 

* 

• 


LO 









16 h ^7 m 24 3 18 5 12 5 D6 3 47 m 00 3 

RA 



Fig. 12. Plot of the positions of the photometrically se- 
lected cluster members (solid circles) and field stars (dots). 
The position of the gap between the two CCDs is indicated 
by the dashed lines. 



OB supergiants either lack photometry due to blending 
in the crowded 'core' regions or are fainter than V^18 - 
the faintest being W15, with V~19 mag. Between V=18- 
19 mag. we have spectroscopy of only 10 of 40 potential 
OB supergiants, although we note that an accurate deriva- 
tion of a limiting magnitude for supergiant candidates will 
necessarily be a function of the significant differential red- 
dening across Wd 1. 

Moreover, we expect that a proportion of the stars 
found between V= 18-20 mag will likely be post-MS but 
of lower luminosity classes - lb or II - and hence intrin- 
sically fainter than the la stars identified to date. Given 
the limitations of the current spectroscopic data set, we 
cannot yet test this hypothesis, but we suggest that the 
the range in V band magnitudes observed for OB stars 
with apparently identical colours (and hence extinctions) 
supports a range of instrinic luminosities for these stars. 

More extreme results are found for the WRs. Only two 
stars are found to have V<18 mag, six from a total of 159 
cluster members have V>18mag. and six are photometric 
non-detections. Moreover, with the exception of the WNE 
star W72, all are either WNL or WCL - if present WNE 
or WCE stars might be expected to be >2 mag. fainter in 
the optical - well below our current detection threshold. 

We therefore conclude that the current spectroscopi- 
cally identified post-MS population for Wd 1, while com- 
plete for the cool transitional objects is significantly in- 
complete for cluster OB supergiants and WRs, ignoring 
the effects of possible binarity and the incompleteness of 
the current photometric dataset due to crowding. Indeed, 
under the assumption of a Main Sequence turnoff at ~07 
V (Sect. 5.1; currently undetectable at V~20.6 mag), a 
large majority of the currently identifiable cluster popu- 
lation will consist solely of post-MS objects. 



Table 3. Comparison of the stellar contents of Wd 1 and 
the three Galactic Centre clusters (Figer 120041 and refs. 
therein). For brevity we define the 'Early Transitional' 
category to encompass sgB[e], LBV and for Wd 1 the 
highly luminous mid-B supergiants W7, 33 and 42. The 
stellar content presented for Wd 1 is restricted to those 
objects for which we have spectral classifications; hence 
while we are likely complete for the cool transitional ob- 
jects the photometric data suggest we are severely incom- 
plete for both OB supergiants and WRs. We have explictly 
excluded the WNL/B Ia+ star W5 and the WNE star 
possibly present as a blend in the W14a spectrum from 
the count of WR stars in Wd 1. Note that the WC con- 
tent of the Quintuplet includes the five Quintuplet Proper 
Members, while for the Galactic Center cluster Genzel 
et al. (j2D03) propose 6 stars with Of/LBV classifications 
which we include as early transitional objects - hence the 
apparent absence of any OBIa stars for this cluster. 

OBIa Earty YHG RSG WN WC~ 
Trans. 

Wd 1 >25 >5 6 3 >7 >6~ 

Quintuplet 14 2 1 5 11 

Arches 20 >6 

Center - 6 2 >10 >10 



Table 4. Comparison of inferred properties of Wd 1 to 
other massive young clusters (after Table 5 of Figer et al. 
1999}. Ml is the mass in observed stars and M2 is the 
total mass assuming a Salpeter IMF between 1-120M Q . 





log Ml 


log Ml 


Radius 


Age 


Cluster 


(M ) 


(M ) 


(pc) 


(Myr) 


Wd 1 


3.8 


4.75 


0.6 


3.5-5 


Quintuplet 


3.0 


3.8 


1.0 


3-6 


Arches 


3.7 


4.3 


0.19 


2-3 


Center 


3.0 


4.0 


0.23 


3-7 


NGC 3603 


3.1 


3.7 


0.23 


2.5 


R136 


3.4 


4.5 


1.6 


< 1 - 2 



5.1. Progenitor masses for the cluster members 

Unfortunately, the current data do not permit the accu- 
rate construction of an HR diagram for Wd 1. This is 
due to the difficulty in accurate spectral classification of 
the early OB supergiants from our present spectroscopic 
dataset. Currently we are unable to constrain the spec- 
tral types of the stars comprising this population to bet- 
ter than 09.5 to B0. 5 la - or a range in temperature of 
6000K, corresponding to an uncertainty in bolometric cor- 
rection of >0.6 mag. A similar problem is also present for 
the population of WRs and is further compounded by the 
lack of photometry for the majority of these stars. We 
further suspect that many of the WRs with detections are 
likely to be binary systems, a problem - compounded by 
blending - that also likely afflicts our sample of OB super- 
giants. Consequently, inferring stellar masses and a cluster 
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Fig. 13. Plot of the mass and radius of Wd 1, selected 
young massive clusters in the Galaxy and LMC and SSCs 
for which dynamical estimates of their mass are available. 
Data from Figer et al. (1999; filled circles), Slesnick et 
al. filled triangles), Smith & Gallagher lj2UUTl open 

triangles), McCrady et al. (20Q2| open stars) and Larsen 
et al. I|2004l open circles). 



age from the current dataset is difficult, particularly due 
to the lack of a detectable Main Sequence. Nevertheless, 
assuming co-evality - justified in Sect. 5.2 - we may make 
progress in determining likely masses for many of the clus- 
ter objects. 

Trivially, the population of WRs are the initially most 
massive objects currently identified in Wd 1. From a study 
of 12 galactic clusters Massey et al. (20.01) report that 
both WNL and WC stars are found in clusters with the 
highest turnoff mass (>50 M Q ), although only one WCL 
star is reported in their study. By contrast, WNE stars are 
identified in clusters with a wide range of turnoff masses, 
with a lower limit of only ^20M©. Following the results of 
Massey et al. (20T!T|) the currently identified WR popula- 
tion of Wd 1 - predominantly consisting of WNL & WCL 
stars - are likely descended from similarly high mass pro- 
genitors. 

Of the transitional stars, comparison of the mean 
luminosity inferred for the YHGs (log(L/L )^5.7) to 
theoretical predictions (e.g. Meynet & Maeder )2003) 
suggests progenitor masses of ^40M©. This is consistent 
with the observation that galactic YHGs occupy a rather 
narrow range of luminosities - log(L/Lo)^5.6-5.8 (e.g. 
Smith et al. 200J) - and hence progenitor masses. In 
this work Smith et al. suggest a possible evolutionary 
sequence for such stars of: 

O MS^BSG^RSG^YHG^WN9-11^WR 

If this scenario is correct, we might also expect the re- 
maining cool and hot transitional objects to share a similar 
mass and luminosity with the YHGs. 



One possible objection to this hypothesis might be 
the apparent presence of only a single LBV in Wd 1, 
a phase which is known to occur for stars in excess of 
log(L/L©)~5.4. However, the harsh UV rich environment 
of Wd 1 is likely to be inimical to the long term sur- 
vival of LBV ejecta, such that only long term monitoring 
may reveal further examples. Obvious candidates include 
the extreme B supergiants W7, 33 & 42a, the WNL stars 
W5 & 44 and BSGs W23a & 43a, which apparently show 
evidence for Ha variability between the 2 epochs of our 
observations. 

Moreover, Smith et al. ( 2004) identify an apparent lack 
of known LBVs in the luminosity range occupied by the 
YHGs and suggest a novel explanation for this - that the 
effect of the bi-stability jump on the mass loss rate of post 
RSG stars evolving bluewards into the S Doradus insta- 
bility strip is sufficient to generate a pseudo-photosphere, 
effectively suppressing the LBV phase for stars of ^40M Q . 
Thus, the apparent lack of LBVs in this luminosity range 
may be a real feature of the evolutionary sequence for stars 
of this mass. 

Finally, as described in Sect. 3.3.1 the presence of 
Ha infilling or emission in the population of OB super- 
giants suggests that they have a mass at the upper range 
- >30M Q - of those observed by McErlean et al. I|1999|l . 
Thus we arrive at a picture of the current observed pop- 
ulation of stars within Wd 1 having progenitor masses 
ranging from 30-40M© (the OB supergiants) through to 
>40M© (the WRs), with the various cool (YHGs & RSGs) 
and hot (LBV and extreme B supergiants) transitional ob- 
jects defining a narrow range around ^40M©. Assuming 
the OB supergiants represent the least evolved members 
of Wd 1 currently identified, a Main Sequence turnoff is 
likely to be found around ~30-35M Q , implying an ^07 V 
classification. As previously described, such a population 
would be below our current detection limit, thus consis- 
tent with the lack of any Main Sequence objects in our 
current photometric or spectroscopic data. 

5.2. The age of Wd 1 

While we may not determine an age for Wd 1 on the basis 
of ioschrone fitting, we note that the presence of the WCL 
population suggests a minimum age of ~3.5 Myrs, while 
the presence of O supergiants implies a maximum age of 
~5 Myr (Meynet & Maeder l2DU5jl . The high luminosity 
and hence mass inferred for the YHGs suggest an age of 
—4 Myrs (Clark & Negueruela 120021 and refs. therein), 
while following Figer et al. ( 2002 ) the presence of relatively 
cool emission line B supergiants implies an age >4Myr. 
We therefore suggest a likely age for Wd 1 in the range 
3.5-5Myr, which is also consistent with our predictions of 
a MS turnoff around 07 V or later. 

An age in excess of ~3 Myr is supported by the 
lack of diffuse radio emission coincident with Wd 1 (e.g. 
Dougherty et al. I2()()4|l . It is expected that for clusters of 
this age a cluster superwind, driven both by the present 
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stellar winds and also a possible supernova component, 
will have dispersed their natal envelopes. Examples of this 
phenomenon include the clusters Danks 1& 2 found within 
the G305 star forming complex (Clark & Porter 2004) 
and HD 32228 found in the LMC Giant Hn region Nil 
(Walborn et al. lT^jl 6 . 

It is difficult to imagine a scenario in which Wd 1 is 
not coeval, as is strongly suggested from the currently 
classified stellar population. The presence of just a few 
very massive stars would be expected to expel the gas not 
yet used in star formation from a cluster and halt fur- 
ther star formation on approximately a crossing time (e.g. 
Goodwin I1997[l . The lack of a significant gaseous compo- 
nent in Wd 1 strongly suggests that this has already hap- 
pened. In order for Wd 1 to be significantly non-coeval 
a mechanism would have to be found that formed only 
low-mass stars before forming all of the high-mass stars 
in a burst. Even if this were the case, our assumption of 
coevality for the observed high-mass stars would still be 
be correct. 

5.3. The mass of Wd 1 

Currently, we have spectroscopically identified ~ 50 stars 
with progenitor masses > 3OM within Wd 1, placing an 
absolute minimum mass ofWdl at 1.5xlO 3 M . Under the 
assumption that all currently identified cluster members 
have masses of >3OM we derive a total mass of observed 
stars some ~4 times larger. Following Table 5 of Figer 
et al. 0999), in Column 2 of Table 4 we compare the 
total mass of observed stars in Wd 1 to those of other 
young compact clusters within the Local Group. With the 
exception of the Arches and R136 clusters, we find that 
Wd 1 is significantly more massive. 

However, such a naked comparison takes no account 
of the either the depth of the relative observations, nor 
the effect that the different ages of the Arches and Wdl 
play. Specifically, at an age of only ^2.5Myrs the mass- 
luminosity relationship is monotonic for the Arches, re- 
sulting in the mass estimate Ml being complete for stars 
>20 M (e.g. Serabyn et al. 1998). In contrast we have 
shown that our current observations are only complete 
for stars with initial masses >3OM in Wd 1. Moreover, 
at an age of ~4Myrs, the mass-luminosity relationship is 
bimodal, such that stars in excess of 4OM will either 

6 We note that a second generation of trigged star forma- 
tion is associated with both the G305 and Nil complexes. 
Examination of both radio (Dougherty et al. 12004ft and mid- 
IR (Midcourse Space Experiment 8-25/^m) data reveals no such 
star formation activity associated with Wd 1; hence we may in- 
fer a lack of significant remnant natal material associated with 
Wd 1. Motivated by the comments of the referee, we speculate 
that either Wd 1 formed from the collapse of the entire na- 
tal molecular cloud (possibly implying a higher than average 
star formation efficiency) or that the superwind generated by 
the unprecedented population of massive stars within Wd 1 
completely dispersed the cloud remnants before a second gen- 
eration of stars could form. 



have been lost to supernovae or will have evolved to the 
WR phase. Indeed, van der Hucht H2001fl demonstrates 
that with the exception of the latest WNLs, WRs of both 
flavours have absolute visual magnitudes M„ > — 5mag. 
- so undetectable in our current data 7 . Thus, we believe 
that our current observations of Wd 1 are only sensitive to 
stars in a narrow range of progenitor masses - specifically 
30-40Mq - and thus may not be directly comparable to 
those of the young <3 Myr clusters such as NGC 3603, 
the Arches and R136 presented in Table 4. 

Following Figer et al. (1993J we present the total mass 
for each cluster in Table 4 assuming an identical Salpeter 
IMF, finding Wd 1 to be a factor of ~2 larger than the 
next most massive cluster, R136. We note however, that 
the current observational position on the full IMFs of mas- 
sive star clusters is unclear. Sirianni et al. (|2000) find that 
R136 has a Salpeter IMF that turns-over at ~ 2M , as 
compared to O.5M0 in the field and local star forming re- 
gions (see Kroupa 2002jfor a review of the IMF). However, 
this has been disputed by Zinneker et al. (2004) who claim 
a Salpeter slope to at least 1M Q . 

Moreover, Figer (200U and references therein) sug- 
gest a significant deviation from a Salpeter IMF for the 
Arches, prompting a downwards revison of the total mass 
to ~ 1O 4 M . The high gas temperatures (even in dense 
molecular clouds) due to the strong ambient radiation field 
in the Galactic Centre may well be responsible for this 
difference from the 'standard' IMF and so the Galactic 
Centre clusters may not be directly comparable to Wd 1 or 
R136 which have formed in more 'normal' environments. 
Indeed, Larscn (2004) finds that the mass-to- light ratios of 
SSCs are consistent with a Kroupa (2002) IMF. Another 
caveat is that these mass estimates are generally based 
on the assumption of virial equilibrium which may be sig- 
nificantly wrong when dealing with very young clusters 
(Goodwin et al. in preparation). 

In summary, it is currently unclear what the form of 
the IMF of massive star clusters is, or indeed if it is uni- 
versal. Sirianni et al. ( 2000 ) do find that the IMF in fields 
close to R136 is Salpeter to at least ~ 1M . It may be 
that primordial and/or dynamical mass segregation acts 
to make the IMF dependent on the size of the region un- 
der consideration, and that abnormal IMFs are the result 
of only considering the cores of clusters. 

Nevertheless, assuming that Wd 1 has a Kroupa two- 
part IMF and that >140 of the current candidate mem- 
bers have masses > 3OM - conservatively defined as the 
number of stars with V<20 mag., excluding the small WR 
component - then we find that its mass must be > 1O 5 M . 
This estimates makes no account of incompleteness and 
should be regarded as a lower limit under the hypothe- 
sis of a standard IMF. Moreover, the additional assump- 
tion that our current observations only sample the stellar 
population within a restricted (~30-40 M ) mass range 

7 We suggest that the photometric detections of a subset of 
the cluster WRs - predominantly the WCLs - may be a result 
of binarity. 
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suggests that we may be significantly underestimating the 
total mass of Wdl. 

A more reliable comparison would be between the 
number count of stars within the 30-40M Q range - cor- 
responding to a mid O Main Sequence star. The Arches 
has ~40 stars in this mass range (Figer et al. 2002; ver- 
sus ~140 for Wd 1. If the IMFs are similar this suggests 
that Wd 1 is at least 3x more massive than the Arches. 
Similarly, Massey & Hunter (|1998|) report that R136 has 
34 stars more massive than 30Mq (My < —4.5) within 
1 pc of the cluster core, and 117 within 10 pc (the nominal 
size of the NGC 2070 cluster). This makes Wd 1 signifi- 
cantly more massive than R136 and at least comparable to 
(and probably more massive than) the NGC 2070 cluster 
at the heart of the 30 Dor complex. It should also be noted 
that the currently known extent of Wd 1 is significantly 
less than 10 pc. 

Similar conclusions may also be drawn for Cyg OB2 
and W49A, both of which have recently been suggested 
to host massive stellar populations. On the basis of near 
IR photometry Knodlseder H2000[> proposed that ~120 O 
stars (or their descendants) are found within Cyg OB2, 
which, depending on the low mass cut off adopted, yields 
a total mass in the range of 4-10 x 10 4 M Q . However subse- 
quent near-IR (Comeron et al. l2002|l and optical (Hanson 
2003 ) spectroscopic surveys were unable to support such 
a conclusion, with Hanson concluding that many of the 
candidate OB stars are likely to be field contaminants 
and hence that the total O star census is unlikely to ex- 
ceed 100. Moreover, at a distance of ~1.7kpc Knodlseder 
(2000) finds Cyg OB2 to have a radius of ^30pc (with a 
half light radius of ^6pc) , an order of magnitude greater 
than we currently infer for Wd 1. 

Similarly, Alves & Homeier |2D03) propose that ~100 
O stars are found within the central 16 x 16pc region of the 
giant Hn region W49A. However, as with Cyg OB2, the 
stellar density appears to be significantly lower than that 
inferred for Wd 1, with only ~30 O stars found within 
the most massive cluster identified (Cluster 1 of Alves & 
Homeier, with a radius of ^3pc). 

Therefore, despite the uncertainties in the appropri- 
ate IMF to apply, Wd 1 appears to be significantly more 
massive and compact than any other local massive star 
forming cluster (Fig. 13). Indeed, the only way in which 
to make Wd 1 comparable to, or less massive than, other 
local massive clusters would be to postulate that it has an 
IMF that turns-over at a very high mass (5-10M Q ) so that 
very little mass is hidden in low-mass stars - arguably an 
equally important result. 

The high mass and extremely compact nature of Wd 1 
make it a clear example of a SSC. Wdl is very proba- 
bly more massive than the average old Galactic Globular 
Cluster (1O 5 M0), with a mass, luminosity and density 
that likely fulfills the criteria expected for proto Globular 
Clusters. 



5.4. Cluster dynamics and evolution 

At least half of the currently observed stellar population 
of Wdl is found within a circular region of radius ~25" , 
corresponding to <0.6 pc at <5.5 kpc. Assuming a mass 
of > 1O 5 M for Wd 1 leads to a minimum density of 
logp=5.OM pc -3 , compared to logp=5.5M pc~ 3 found 
for the Arches cluster. Given that the distance is likely 
to be less than 5.5 kpc and that this estimate does not ac- 
count for incompleteness (which is expected to be worst in 
the centre of the cluster), we consider it likely that Wd 1 
is of comparable density to the Arches. 

The collisional timescale for stars in the centre of Wd 1 
is of order 10 4 Myr (Binney & Tremaine ll987f) assuming 
a central number density of stars of a few x 10 6 pc -3 for 
solar-type stars. This number will be significantly reduced 
by the presence of many large (in both mass and radius) 
stars (see for example the simulations of Porteies-Zwart 
et al. I2004|) . We might expect a collision between two of 
the stars in the core of Wd 1 every 10 3 — 10 4 yrs, and so 
it is not unreasonable to propose that W9 is an unrelaxed 
recent merger remnant. 

Recent work by Portegies-Zwart et al. (j2004[) suggest 
that under such extreme conditions runaway stellar merg- 
ers may lead to the production of an intermediate mass 
black hole. Accretion onto such objects has been suggested 
by some authors to explain the class of Ultraluminous 
X-ray Sources observed in external galaxies. // such an 
object is present in Wd 1, the lack of such a luminous 
(L x > 10 39 erg/s) X-ray detection (Clark et al. I1998jl may 
indicate that either no mass donor is associated with it, 
or that mass transfer is either transitory or proceeds at a 
comparatively low rate (e.g. via direct wind fed accretion 
rather than Roche Lobe overflow). 

Assuming a velocity dispersion in the core of ~ 15 
km s _1 (expected if the cluster core is approximately viri- 
alised) would imply that the core is > 100 crossing times 
old. This would mean that the core is dynamically re- 
laxed. It is possible that the outlying stars visible in Fig. 
12 have been dynamically ejected from the core. Indeed, 
it is quite possible that some stars may have travelled in 
excess of 50pc from the cluster centre. Whilst the core of 
Wd 1 should be dynamically relaxed, the outer regions 
are probably dynamically young and substructure in the 
initial cluster may not have been erased (e.g. Goodwin & 
Whitworth l2004|l as is thought to be the case in the NGC 
2070 cluster fMevlan HS55|l . 

The huge population of massive stars in Wd 1 suggests 
that it will have a very significant impact on the local ISM. 
We estimate that the total number of type-II supernovae 
(SNII) will exceed 1500. This implies that Wd 1 will inject 
> 10 55 ergs of energy from stellar winds, UV radiation and 
SNII explosions on a timescale of < 40 Myr. The energy 
flux from Wd 1 will exceed 10 40 ergs s _1 for several Myr. 
The lack of radio emission and an HII region associated 
with Wd 1 (Sect. 5.2) suggests that Wd 1 already con- 
tains very little gas, presumably due to feedback from the 
massive stellar population. Indeed, the impact Wd 1 has 
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already had upon its environment may be the reason why 
it has escaped identification as such an extreme cluster for 
so long. 

The input of energy from many hundreds of SNII 
within Wd 1 would be expected to drive a significant 
superbubble and evacuate a large portion of the local 
Galactic disc (until destroyed by the differential rotation 
of the Galaxy). Wd 1 could well eject a significant amount 
of material out of the Galactic disc, driving a fountain or 
creating high velocity clouds. 

Wd 1 is also expected to return > 1O 3 M of heavy ele- 
ments to the ISM (e.g. Tsujimoto et al. 11995(1 . Due to the 
extremely high mass of many of the stars in Wd 1 com- 
pared to the content of lower-mass clusters, Wd 1 would 
be expected to return unusual abundances of heavy ele- 
ments (e.g. Goodwin & Pagel l20Ul) l. 

6. Conclusions 

We have presented both spectroscopic and photometric 
observations of the cluster Wd 1. On the basis of their 
optical colours, we have identified ~200 candidate mem- 
bers, of which 95% are found to lie within 2arcmin of the 
nominal cluster core, and ^50% within an inner region 
of radius ~25" . Adopting an upper limit to the distance 
of 5.5kpc - from consideration of both spectroscopy and 
photometry - therefore implies that Wd 1 is rather com- 
pact; the angular radii corresponding to <3.2 and <0.6pc 
respectively. 

We have obtained spectral classifications for 53 candi- 
date cluster members, all of which are found to be post 
MS stars. Of these we find 14 to be WR stars of both 
flavours, 25 to be OB supergiants with spectral types of 
approximately BOIa, and at least 14 transitional objects. 
The latter group consists of both high (LBV, sgB[e] & 
extreme BSGs) and low (YHGs and RSGs) temperature 
objects. In doing this we confirm the result of West87 in 
identifying 4 YHGs within Wdl and add a further 2 ex- 
amples, leading to a total population of such stars that is 
directly comparable to that of the Galaxy or the LMC. 

Unfortunately, we are currently unable to construct an 
HR diagram for Wd 1. This is due to the uncertainties in 
the spectral classification of the OB supergiants, which 
preclude an accurate determination of their intrinsic lu- 
minosities, compounded by significant differential redden- 
ing across the cluster and the presence of a non standard 
extinction law. Nevertheless, we may use spectroscopic 
diagnostics to infer a mean luminosity for the YHGs of 
log(L/L©)^5.7, placing them firmly at the HD limit and 
implying ~40M Q progenitors. Following the recent work 
of Smith et al. H20Q4J we infer similar initial masses for 
the other transitional objects within Wd 1. While no un- 
ambiguous spectroscopic luminosity diagnostics exist for 
OB supergiants, following the results of McErlean et al. 
( 1999$ we infer progenitor masses of ~30-40M Q from the 
presence of Ha emission in their spectra. 

These estimates suggest that we might expect a Main 
Sequence turnoff mass for Wd 1 of around 30-35M©, con- 



sistent with a population of late O dwarfs with spectral 
types no earlier than ~07 V. Adopting a median value 
for V— My ^25.3 mag. implies that such objects will be 
found at V~20.5 mag. This is entirely consistent with the 
lack of MS objects in our spectroscopic dataset, and the 
apparently lack of a MS in the photometric data. This in 
turn implies that a large majority of the ^200 currently 
identified cluster members are also likely to be post-MS 
objects. 

Simply considering spectroscopically classified stars 
within Wd 1, we find a minimum mass of the order of 
1.5xlO 3 M0. However, the twin assumptions that ^140 
currently identified cluster members had initial masses 
>30M Q and that a Kroupa fMtyl two-part IMF is ap- 
propiate for Wd 1 yield a total cluster mass of ^10 5 M Q . 
Note that this estimate takes no account of incomplete- 
ness in the current datasets nor that they are likely only 
sensitive to stars within the restricted ~30-40Mq range, 
both of which will act to increase the mass of the clus- 
ter. Nevertheless, Wd 1 appears significantly more massive 
than any currently identified open cluster in the Galaxy, 
and indeed is likely to be more massive than the average 
Galactic Globular Clusters. We conclude that Wd 1 rep- 
resents the first known example of an SSC in the Local 
Group - i.e. a factor of ^1000 times closer than the next 
closest SSC, NGC1569-A (Hunter et al. 1277071)) . 

With such an extreme mass and population of massive 
stars Wd 1 promises to greatly advance our understanding 
of the formation and evolution of massive stars, both in- 
dividually, and in the extreme physical conditions present 
in SSCs. Given the wealth of post-MS objects, the identi- 
fication of a MS and construction of an HR diagram will 
allow accurate progenitor masses to be assigned to differ- 
ent spectral types, permitting the refinement of current 
theoretical pathways through the post-MS 'zoo' and the 
empirical verification of the HD limit for a population of 
co-eval stars at a single metallicity. Deep adaptive optics 
observation will permit the identification and study of sub- 
solar mass objects within Wd 1 - if their formation has not 
been inhibited by the population of massive stars. 

Indeed the presence of Wd 1 within the galaxy will 
provide us with a unique insight into the physical pro- 
cesses occurring in the SSCs heretofore only identified in 
external star burst galaxies, challenging both our current 
theories of star formation and the role it plays in the wider 
evolution of the Galaxy. For example - assuming a nomi- 
nal star formation efficiency of ~10% - how did ^10 6 M Q 
of gas collect in a region ~one pc across? The average den- 
sity of the pc 3 molecular cloud from which Wdl formed 
must have been in excess of 10 6 atoms cm -3 and hence 
the star formation process in this region must have been 
very different from that in local star forming regions such 
as Taurus or Orion. Recent simulations by Bonncll et al. 
(2004) of clustered star formation in a highly turbulent 
molecular cloud give probably the closest insight into star 
formation in Wd 1. 
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Appendix A: OB supergiant classification criteria 

Caron et al. ( 2003;) have recently outlined a classification 
scheme for 09 to B5 la - V stars between ~8400 to 8900A, 
based on both observations of standard stars and synthetic 
models. We have followed this methodology in the exten- 
sion of this scheme to early O - late B supergiants, via the 
construction of an extended grid of non-LTE synethetic 
spectra, normalised where possible to observations of OB 
supergiants. We have further extended this to include the 
6500-7900A range to encompass our full data set, although 
we concentrate on the I band for our primary classification 
criteria. 

We calculated our grid of OB supergiant models with 
the non-LTE atmospheric code CMFGEN (Hillier & Miller 
HMty between T e// = 10kK and T e// =50kK at 2000K in- 
tervals, accounting for line blanketing by H, He, C, N, 
O, Mg, Al, Si, S, Fe. Our approach follows the approach 
outlined for OB stars by Crowther et al. (J2002), with an 
assumed 20kms~ 1 microturbulence. A uniform rotational 
broadening of 75kms~ 1 is adopted throughout. 

The temperature scale for O supergiants has recently 
been revised downward as a result of the incorporation 
of line blanketing (Crowther et al. 1200 2 Herrero et al. 
2003.; Repolust et al. l2004[> . Similarly, the scale for early- 
B supergiants follows Crowther et al. (2004) whilst mid 
and late-B stars follow recent unpublished work by one of 



us (PAC). Luminosities for individual models follow from 
theoretical bolometric corrections, plus the Humphreys & 
McElroy (1984$ absolute magnitude-spectral type calibra- 
tion. 

We have adopted wind velocities for each subtype fol- 
lowing Prinja et al. I|1990(l and Lamers et al. (|1998f) . and 
mass-loss rates following the wind-momentum relation- 
ships from Kudritzki et al. I|1999[l for B stars 8 and Puis et 
al. l(TM>|l for O stars. 

The synthetic spectra in the rang e 8200 to 8900A are 
plotted in Fig. 13 with, where available, standard spec- 
tra overplotted for comparison (the spectra of standard 
stars are summarised in Table 5). A summary of the 
model parameters for each spectrum, along with optical 
and near IR colours and V band bolometric correction are 
presented in Table 6. Comparison of synthetic spectra to 
observational data indicate an encouragingly close corre- 
spondance, given that models have not been tailored for 
the individual standard stars. 

The first requirement of an I band classification scheme 
for the current dataset is to distinguish between Main 
Sequence and Supergiant spectra, which Caron et al. 
( 2003 ) demonstrate may be accomplished via the FWHM 
of both the Paschen series and Hei absorption features, 
which are systematically narrower for the supergiants 
due to the Stark effect (their Fig. 6). As a consequence, 
Paschen lines shortwards of Pal6 and the weak Hei tran- 
sitions such as 8775A are unresolvable in Main Sequence 
stars, while clearly present in supergiants of similar spec- 
tral type. 

Accurate spectral classification of OB stars is more 
problematic, with Caron et al. (2()()3i forced to rely solely 
on the absolute strengths of the Paschen lines for classifi- 
cation of stars in the range 09-B5. While we too are forced 
to partially adopt this approach, particularly for stars with 
temperatures in excess of 34kK (071), we are able to iden- 
tify a number of further temperature diagnostics in both 
standard and synthetic spectra. In attempting to describe 
a classification scheme we adopt the philosophy of rely- 
ing on observational over synthetic data, and the occur- 
rence of specific lines and particular line ratios rather than 
the absolute line strengths. For a line to be identifiable 
we require an EW>0.2A in either our synthetic or stan- 
dard spectra. In practice all lines listed satisfy both cri- 
teria, with the exception of the Pal6+Cni~8500A blend, 
which is absent in our 081a (34kK) synthetic spectrum 
but clearly present in the standard spectra. Hen 8238 A 
lies outside the range of the standard spectra presented, 
but its occurrence is confirmed in unpublished O star spec- 
tra of one of us (Crowther). 

Detailed criteria for spectral types 02-AO are: 
•02 to 07 la (50 to 34kK): The spectra of the earli- 
est stars are dominated by the presence of He n 8238A and 
Paschen series absorption lines alone. Consequently accu- 

8 We adopted an identical form of the wind momentum rela- 
tionship for late B supergiants to that established by Kudritzki 
et al. •: l!)!lTJJ for mid-B supergiants. 
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rate spectral typing within this temperature range must 
rely solely on the relative strengths of the Paschen series, 
which our synthetic spectra indicate increases by a factor 
of ^2 between 02 to 07. Without adequate observations 
of standard stars to normalise our synthetic spectra to, 
we refrain from providing an absolute calibration of line 
strengths against spectral types (noting that fortuitously 
this does not affect our conclusions for Wd 1). 

•08 to 09 la (34 to 32kK): While superficially 
similar, spectral types 08 to 09 may be distinguished from 
hotter O stars by the emergence of weak He I absorption 
features, most notably at 8775A. An additional diagnostic 
is the Pal6+Cin blend at ^8500A, which has a greater 
depth and strength than the adjacent Pal5 feature. 

•09.5 to BO la (30 to 28kK): 09.5 to BO la spec- 
tra are characterised by the abrupt disappearance of He II 
absorption. As with 08 to 09 la stars, the Pal6+Cin 
~8500A blend exceeds the strength of Pal5. 

•BO. 5 to B1.5 la (26 to 20kK): Spectra consist 
solely of Hi and Hei features; as with 02 to 07 la stars, 
detailed sub classification relies on the absolute strength 
of the Paschen series (c.f. Caron et al.|2003jfor a detailed 
analysis). 

•B2 (18kK): The appearance of the Oi 8446A line 
with a central intensity less than the adjacent Pal 8 line 
acts as a diagnostic for B2 la stars. 

•B3 (16kK): As B2 la, but marked by the emergence 
of weak Ni absorption features in the range ~8670 to 
8730A. 

•B5 (14kK): As above, but the intensity of Oi 8446A 
now exceeds that of Pal8. 

•B8 to AO (12 to lOkK): As B5 la but with the loss 
of residual Hei absorption features, leaving a spectrum 
dominated by Paschen, Oi and Ni absorption. 

Finally, an additional classification criterion found in 
the 6500 to 7900A band is the O I 7774A blend, which 
is detected (EW>0.6A) in both standard and synthetic 
spectra of spectral types of BO. 7 and later. Therefore, if 
available, this additional diagnostic helps to break the de- 
generacy found for the 26 to 20kK (BO. 5 to B1.5Ia) tem- 
perature range. 



Table A.l. Individual standard stars used for the con- 
struction of Figs. 4-5. Data are from 1 Munari & Tomasella 
(QUEl, 2 Le Borgne et al. (SHU, and 3 Cenarro et al. 
(2001J. All stars are either type la or lb supergiants, the 
latt er given in itallics. 



Spectral type Star ID 

04 HD 190429 1 

06 HD 210839 1 

07 HD 192639 1 , HD 57060 3 

08 HD 167971 1 

09 HD 210809 1 , HD 5706 1 3 

09.5 HD 228779 1 ,HD 306U\HD 37742 1 , HD 209975 3 

BO HD 37128 1 

B0.5 HD 194839 3 

B2 HD 268623 2 , HD 206165 3 
B3 HD 198478 1 , HD 271163 2 , HD 14134 3 

B5 HD 164353 1 ' 2 , HD 13267 3 

B8 HD 34085 1 , HD 199478 3 

AO HD 87737 1 , HD 39970 3 
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Fig. A.l. I band (8200 to 9000 A) synthetic spectra between 10 and 50 kK (solid black lines; parameters given in Table 
3). Prominent transitions of Hi, Hei, Hen, Oi, Ni and Can indicated. Published OB supergiant spectra overplotted; 
data from Munari & Tomasella l|1999l red dotted lines), Cenarro et al. (20(12 blue dotted lines) and Le Borgne et al. 
(2003 cyan dotted lines); stellar identification given in Table 2. 

Table A. 2. Stellar properties adopted for the construction of the grid of synthetic OB supergiant spectra presented 
in Fig. 3, presented with the resultant absolute V band magnitude, (V — I) and (V — K) colour indices, and V band 
bolometric correction (B.C.). 





Spec 


log(L/L ) 


logM 




M v 


(V-I) 


(V-K) 


B.C. 


(kK) 


Type 




(M@ yr" 1 ) 


(kms -1 ) 


(mag.) 


(mag.) 


(mag.) 


(mag.) 


50 


02 


6.24 


-4.80 


3000 


-6.5 


-0.41 


-0.84 


-4.35 


46 


03 


6.13 


-4.92 


2600 


-6.5 


-0.41 


-0.85 


-4.07 


40 


04 


5.95 


-5.18 


2400 


-6.5 


-0.39 


-0.83 


-3.64 


38 


05 


6.02 


-5.09 


2300 


-6.8 


-0.40 


-0.81 


-3.49 


36 


06 


5.94 


-5.20 


2200 


-6.8 


-0.40 


-0.87 


-3.31 


34 


08 


6.04 


-5.08 


2100 


-7.2 


-0.40 


-0.88 


-3.15 


32 


09 


5.97 


-5.18 


2000 


-7.2 


-0.41 


-0.89 


-2.97 


30 


09.5 


5.89 


-5.26 


1900 


-7.1 


-0.42 


-0.88 


-2.87 


28 


BO 


5.75 


-5.8 


1750 


-7.0 


-0.38 


-0.84 


-2.62 


26 


B0.5 


5.69 


-5.9 


1600 


-7.0 


-0.37 


-0.82 


-2.47 


24 


B0.5 


5.60 


-5.9 


1200 


-7.0 


-0.35 


-0.82 


-2.25 


22 


B0.7 


5.52 


-6.0 


1000 


-7.0 


-0.33 


-0.74 


-2.05 


20 


B1.5 


5.43 


-6.9 


800 


-7.0 


-0.32 


-0.70 


-1.83 


18 


B2 


5.49 


-6.7 


550 


-7.4 


-0.29 


-0.63 


-1.59 


16 


B3 


5.30 


-6.9 


400 


-7.2 


-0.22 


-0.55 


-1.30 


14 


B5 


5.09 


-7.2 


300 


-7.0 


-0.14 


-0.46 


-0.97 


12 


B8 


4.93 


-7.4 


250 


-7.0 


-0.05 


-0.34 


-0.57 


10 


AO 


4.76 


-7.3 


200 


-6.73 


+0.07 


-0.17 


-0.15 



